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ABSTRACT 

The sheath current -voltage characteristics for a plane cathode emitting current from 
a thermionic arc spot are derived. The derivation is based upon the requirement of en- 
ergy conservation at the cathode spot and also in an ionization region adjacent to the spot. 
The results for given gas conditions and electrode work function predict a sheath voltage 
which decreases with both increasing current for a fixed undisturbed (i. e. , far from the 
spot) cathode temperature, and increasing undisturbed cathode temperature at a fixed 
current. Experimental data is shown which supports the conclusions predicted by the 
model. 
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ANALYTICAL AND EXPERIMENTAL STUDIES OF MHD GENERATOR 
CATHODES EMITTING IN A "SPOT” MODE 
by Lester D. Nichols and Maris A. Mantenieks 
Lewis Research Center 

SUMMARY 

A simple physical model for the behavior of thermionic arc spots on MHD generator 
cathodes is presented. In this model ions are accelerated through a sheath to the cathode 
by the cathode voltage drop. Both kinetic energy and ionization energy are transferred to 
the cathode. A spot forms whose size and temperature are determined from an energy 
balance between the heat transferred to the cathode by the ions (and possibly excited state 
atoms) and the heat transferred from the spot by conduction in the cathode. 

It is assumed that the ions are created by electron collisions in the gas in a region 
near the cathode. The electrons are thermionically emitted from the cathode. They are 
accelerated from the cathode into this ionization region by the same cathode voltage drop 
which accelerates the ions to the cathode. An energy and mass balance in this ionization 
region along with the equation of state for a perfect gas can be used to determine the ion 
current. 

Based upon the present model, the effects of varying the current, gas pressure, and 
undisturbed (i. e. , far from the spot) cathode temperature are predicted. An increase in 
any one of these parameters (while holding the other two parameters constant) will result 
in a decrease in the cathode sheath voltage drop and an increase in the cathode spot tem- 
perature. The cathode spot size and heat transfer to the cathode spot increase with in- 
creasing current, but decrease with increasing pressure and undisturbed cathode temper- 
ature. This model can be used in MHD generator design to predict the effects of varying 
load current, gas pressure, and electrode temperature upon cathode spot parameters. 

The electrical characteristics of 2 percent thoriated tungsten cathodes used in a 
supersonic MHD generator have been studied experimentally. The generator working 
fluid is pure argon which is arc heated. A single pair of electrodes is studied. The elec- 
trodes are flat plates held parallel to the flow direction and away from the generator wall. 
Illuminated spots are apparent on the cathode when current is flowing in the generator 
load circuit. At the same time a region of higher ion concentration in the argon near the 
spot is observed by spectroscopic measurements. The cathode voltage drop and temper- 
ature of the spot have been measured as a function of the arc spot current. The results 
of the experiment are in agreement with the predictions of the model. 



INTRODUCTION 


The mode of electron emission from an MHD generator cathode can affect the per- 
formance of the generator. Under certain conditions the emission of the electrons is 
confined to a small spot on the cathode. These spots are similar to the spots observed 
in any thermionic arc (ref. 1). There are many theoretical studies of these arc spots, 
and they have in common the assumption that the spot is formed by ion energy transfer to 
the cathode. The ions are formed in a region near the cathode as a result of electron 
collisions with gas atoms. The differences among the previous studies involve their 
treatment of the ionization region and the number of unknown parameters which must be 
specified in order to use the model to predict current voltage characteristics. Ecker 
(ref. 2) determines the ion current from the solution which gives the minimum voltage as 
a function of a ’'contraction parameter" rather than considering the energy conducted 
from the spot. Bauer (ref. 3) and Lee and Greenwood (ref. 4) leave the ion current as 
a free parameter to be determined experimentally. Adams and Robinson (ref. 5) con- 
sider a rather complete set of heat transfer mechanisms, but require the specification of 
the angle between the current stream tube and the normal to the cathode before compari- 
son between their model and experiment can be made. The aim of this study is to pre- 
sent a model which requires no empirical constants. Such a theory is presented by 
Bade and Yos (ref. 6). But they make many assumptions about the nature of the phenom- 
ena to develop a simple theory. The present report is an attempt to improve the Bade 
and Yos theory by using fewer assumptions but retaining simplicity. Bade and Yos as- 
sume that the ions and electrons are at the same temperature, all the atoms are ionized, 
and that the electron current carries no energy from the ionization region into the dis- 
charge. We offer a means of calculating electron temperature, ion density, and the en- 
ergy carried by the electrons from the ionization region to the discharge volume. In all 
other respects these two models are the same. 


THEORETICAL MODEL 
Cathode Energy Conservation 

The electron current in a thermionic arc (i. e. , an arc whose electron current can 
be accounted for by thermionic emission from the cathode) is supplied by thermionic 
emission from a spot of radius R on the cathode. The size of this spot can be estimated 
from an energy balance at the cathode surface. The energy impinging on the surface is 
assumed to come from an ion flux. This ion flux, in turn, is assumed to originate in an 
"ionization region" shown in figure 1. The ionization is accomplished as a result of the 
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Figure 1. - Sketch showing cathode spot, space charge region, ionization region, 
and discharge region. Also shown are descriptions of type of assumed 
phenomena occurring in each region. 

acceleration of the electrons emitted from the spot. The energy transported by the ion 
flux is proportional to the ion current times the net potential energy of the ions (ref. 1): 

i i {y sh + ( Pi-<Pj 


This expression assumes neutralization of all ions at the cathode (no ion reflection), and 
all of the ion energy is transferred to the cathode (an accommodation coefficient of 1). 
There is also a flux of excited neutrals to the cathode. The energy transported by this 
flux is j This expression assumes deexcitation at the cathode. The electron cur- 
rent emitted from the spot takes with it energy: 


Herein we will consider work functions such that the electrode spot temperature contribu- 
tion can be neglected. 

If the spot area is small compared to the electrode area, then the energy conducted 
from the spot is (ref. 9) 
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4RK(T s - T J 


If the cathode is in radiation equilibrium with its surroundings at zero current (i. e. , with 
no spot), then the additional energy required to maintain the additional radiated energy 
because of the spot is 


7rR 2 ea(T^ - T*) 

If there is no mass transfer from the cathode, then the steady state energy balance must 
be 

4RK(T S - T J = »R 2 [ji(V sh + n - <e w > + ll^l ' Vw - £<r ( T s - T i)] <l) 

This equation provides a relation between the spot temperature, spot radius, and the 
various particle currents. The electron current can be related to the spot temperature 
by the Richardson equation 


•*e ^rich^s 


exp - 



( 2 ) 


The ion current can be related to the electron current if one considers a mass balance 
and equation of state in the ionization region. 


Ionization Region 

Mass conservation . - Assume that the ion current is given by the random thermal 
current: 


where 


ii = 4 eN i v i 



(3) 


( 4 ) 
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The ion temperature in the ionization region is assumed to be equal to the neutral temper- 
ature, and both are assumed to be the same in the ionization region as in the discharge 
region. The ion number density will, of course, be influenced by the electrons which are 
carried to the region, hi order to calculate this quantity, an ionization model must be 
established. We will assume that the atoms have three energy levels: ground state, first 
excited state, and ionization. Also, we assume that the pressure in the ionization region 
is the same as in the remainder of the gas. This assumption will be accurate only when 
the pressure gradients required to provide the particle flux lead to small pressure 
changes in the ionization region. For this case, then 


— s N = N + N 1 + N- + — N 
kT ° 1 1 T e 


( 5 ) 


since the atoms and ions all have the same temperature. If one also assumes charge 
neutrality and defines the temperature ratio r as 


rs ^e 


( 6 ) 


and the parameter p as 


, N o + N l 
N i 


(?) 


then equation (5) can be written as 


N i 


N 1 + r + p 


( 8 ) 


To estimate the size of the parameter p, consider the continuity equations which deter- 
mine the number density of the three levels of the model atom. A three level model 
atom oversimplifies the actual situation, but the results can be used to determine trends 
(ref. 7). The collision cross section for excitation from the ground state is Q q j, and 
the collision cross section for ionization from the excited state is Q-q and from the 
ground state is Q q .. The continuity equations can be integrated over a planar ionization 
region of thickness X (see fig. 1) and written as 
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ii = Xe (<Q 0i v e N e >N 0 + <Q u v e N e > Nl ) 

il = Xe (<Q ol v e N e >N D - <Q u v e N e )H 1 ) 

>o - -Xe (<Q oi v e N e > + <Q ol v e N e >)N 0 

The variables have been assumed to be constant in the ionization region and ( ) denotes 
an average over the electron energy distribution. These continuity equations imply no 
radiative deexcitation or recombination in the ionization region. This is a good assump- 
tion for the neutral pressures and electron temperatures considered in this study (refs. 7 
and 10). 

The probability of ionization from the excited state is larger than for the other two 
collisions considered, assuming that once an atom is excited to this state it will be ion- 
ized immediately, and that the flux of excited particles from the ionization region will be 
zero. This assumption is valid for a large ratio of electron current to flux of ground 
state atoms. Then, equation (10) implies 


(9) 

( 10 ) 
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Equation (9) can now be written as 


h = XeN o ( <Q oi v e N e > + < Q ol v e N e>) 

If one substitutes for j. from equation (3), equation (13) becomes 

\ v i N i = XN o (<Qol v e N e > + H^eV) 


( 12 ) 


(13) 


(14) 


The thickness of the ionization region X must be determined. Herein, we assume 
that this thickness is equal to the electron mean free path for collisions which result in 
the excitation of a ground state atom to the first excited state. The mean free path for 
ionization of these excited-state atoms is shorter than that for excitation, so that ioni- 
zation from the excited state can readily occur in this ionization region. The mean free 
path for ground state ionization is larger than for excitation, but for the sheath voltage 
range considered herein (less than twice the ionization potential) this process is less 
likely than excitation. At higher voltages, however, one should use a different length to 
determine the ionization region thickness. Thus, in this study, we take X to be 
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(15) 


With this definition of X and using equations (12) and (14), equation (7) may be written as 


- v.N. 


i o 


P = 


<«li v e N e>^5olVy ] / <Qol v e M e> 
tvO \<Qoi v e N e> + <«*ol v e N e>/ \ <«li v e N e> 

G G 


(16) 


In order to evaluate the term p, the electron distribution function must be known. It 
is quite likely that the electrons will be maxwellized (ref. 8) at some temperature less 
than the sheath voltage. We have used Dugan and Sovie results for electron temperatures 
from 1 to 25 volts and have found that the product of the terms multiplying the current 
ratio in equation (16) is less than 1. Therefore, if the electron current is large com- 
pared to the thermal flux of ground state atoms and if the temperature ratio r is large 
compared to 1, then p can be neglected. This implies that the degree of ionization in 
the ionization region is greater than 30 percent. Assuming this to be the case, equation 
(13) becomes 


N i 


N 


1 

1 + r 


(17) 


Next, the electron temperature in the ionization region must be determined. We as- 
sume that the average energy of the electrons is less than the energy required for ioni- 
zation from the first excited state, since this is the inelastic collision with the smallest 
energy transfer from the electron to the atom and, hence, the smallest mean free path. 
One may assume that the distribution of electrons with energy less than this energy 
(<p. - cp^) will be such that the average energy will be one-half of this energy. That is, 
the electron temperature in the ionization region is 

-kT e = -e(<p. - cpA (18) 

2 e 2 1 1 

Therefore, the ratio of temperatures can be written as 

r _ _ e <n - 

T 3kT 
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Bade and Yos perforin the calculations in their model with r = 1. The above choice of 
r is the main difference between the present model and that of Bade and Yos. 

Energy conservation . - An energy balance in the ionization region is 

je V sh = h*i + (j e + W (19 > 

where all of the ionization energy is assumed to return to the cathode, and the electrons 
leaving the ionization region carry the average energy into the discharge volume. The 
ratio of ion to electron current can be written as 



v sh - \ - n> 

•Pi + \ <<°i - '"l) 


( 20 ) 


The cathode spot energy equation (1) can be written using equation (19) as 


4K(T g - T J = ttR^O 


jt>i + ^[Vsh " ^w 4 ( ^ ' " £a ( T s ' T i) 


( 21 ) 


If radiation is negligible (i. e. , when the term involving the fourth power of temper- 
ature is small compared to the other terms in brackets), one can introduce 


1 = tfR 2 (i e + jj) 


( 22 ) 


and rewrite equation (21) as 


4K(T S - T J = * XT 
V ff(j e + 

Introduce the Richardson equation (eq. 


4K(T g - T J = flA rlch exp 


7j °e + ii>| v sh ’’l*} 


(2)) and from equation (20) 



(23) 


(24) 
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A function $ can be defined which does not depend upon the electrode temperature as 


$ = + 1 + 


V sh “ 2 (<P i " ^ j(V sh - cp w ) j ((Pi - <?!> 


n + \ ^i " *1> 






(25) 
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(26) 


Equation (24) can then be written as 




® i i 

± exp — 

I ref * 0 20 


(27) 


This gives the current as a function of the spot temperature and the sheath voltage. 
Another relation between sheath voltage and temperature is determined by a combination 
of equations (3), (17), and (20): 


v S h - - n> A T 2 cxp 

<?i + -(<?i - <P X ) riCh S ^s-l 41 + ] 
2 


(28) 


The ion current can be written in terms of the gas pressure and temperature as 


9 



•*i, max 


eP 


(29) 


eNVj 

4 


i 


2kT jrm. 


;pVr 


V 2kT e 


Equation (28) can then be written as 



V sh - - fl') 


1 + r 

V~r 



(30) 


Define a reference pressure as 


ref 


A rich 




12 2 

For argon and a cathode work function of 3. 0 eV, P re f = 2. 824x10 N/m . Then, if one 
defines a function as 


* =. 


1 + r 

P V sh " - (< ?i - ?1> 

F ref sn 2 1 1 

/ 

P Vi + kVi-Vt) 


(31) 


equation (30) becomes 


= - exp — 

e 2e 


(32) 


This provides a relation between the gas properties, the work function, and the spot tem- 
perature. The function is plotted in figure 2. 

To use this model one needs to specify only the gas, its pressure and temperature, 
the electrode work function, and undisturbed cathode temperature (away from the spot) in 
order to calculate the current -voltage characteristics. The radius of the spot, the heat 
transferred to the spot, and the current density can also be calculated. 

The radius of the spot is 


R = 



_I 

e + 


(33) 
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Dimensionless cathode spot temperature, 0 

Figure 2. - Function of cathode sheath voltage and ion current density detmed 
in equation (32) plotted as function of dimensionless cathode spot temperature 


This can be rewritten using equations (2), (20), and (26) to read 



where 



The heat transferred to the spot is 


Q = 4RK(T g - TJ 


(36) 


which can be rewritten with the aid of equations (27) and (34) to read 

Q _ J I \ $ 

^ref \ ref/ |/l + ^ 


(37) 


where 


ecp 

«ref ’ R ref K 


(38) 


The current density can be determined from equations (2) and (20) 


J = J e + h = A rich T s 


exp - 


/ecp. 


w 


, kT, 


(1 + rj) 


(39) 


which can be rewritten using equation (32) to read 

= (l + yp — (40) 

^ ref ^ 

where 

Vef * \ich (~) <«> 

These parameters, along with the sheath voltage and spot temperature from equa- 
tions (27) and (32) are shown in figures 3 to 5 for a work function of 3. 0 eV. In each fig- 
ure these parameters are shown for two ion temperatures, 500 and 1000 K, and for 
12 

P ref = 2 - 824x10 newtons per square meter. In figure 3 the current is varied; in fig- 
ure 4, the cathode temperature; and in figure 5, the pressure. As the current increases 
the sheath voltage decreases, and the other four quantities increase as can be seen in 
figure 3. The heat transfer increases rapidly with increasing current. 

The sheath voltage and heat transfer decrease with increasing temperature; the 
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meter; undisturbed cathode temperature, 1800 K. 


Ion temperature, 

T|. 

K 



Undisturbed cathode temperature, T^, K 

Figure 4. - Cathode spot parameters as function of un- 
disturbed cathode temperature. Work function, 3.0 
electron volts; pressure, 3. 5x1$ newtons per square 
meter; current, 100 amperes. 




Figure 5. - Cathode spot parameters as function of pressure. Work function, 3. 0 electron volts,* current, 
100 amperes; undisturbed cathode temperature, 1800 K. 


other three quantities increase as can be seen in figure 4. The heat transfer is rather 
sensitive to the temperature. 

In figure 5 it can be seen that the current density increases very rapidly with in- 
creasing pressure while the spot radius decreases and the spot temperature increases. 
The sheath voltage and heat transfer pass through a maximum value as the pressure in- 
creases. The low pressure region is characterized by a big spot and rather low temper- 
ature. The high pressure region is characterized by a small spot and high temperature. 


EXPERIMENTAL APPARATUS AND MEASUREMENTS 

Experiments conducted in a nonequilibrium argon plasma MHD generator led to the 
development of the physical model to explain the cathode spots. The experiment consisted 
of studying applied or magnetically induced voltages across a pair of plane -parallel elec- 
trodes immersed into hot -flowing argon. The electrodes were immersed to raise their 
temperature (T^ in the preceeding analysis) to provide thermionic emission. 
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Description of Apparatus 

Arc heater and nozzle . - The arc heater is vortex stabilized and consists of a 
li-inch- (32-mm-) diameter, annular, water-cooled copper anode and a water-cooled 

O 

thoriated tungsten cathode. The argon flows into the arc through tangential holes. 

A water-cooled convergent-divergent copper nozzle with a throat diameter of 1. 44 
inches (3. 65 cm) and an area ratio of 4. 34 is used. The Mach number at the nozzle exit, 
calculated for isentropic expansion of a perfect gas, is 3.6. 

Test section. - The test section is a 7. 6 -centimeter -inside -diameter stainless steel 
tube 35. 6 centimeters long. Thoriated tungsten (2 percent ThOg) plate electrodes 
(0. 76 mm x 1 . 27 cm x 2. 22 cm) with pure tungsten stems are placed into the stream. 

To measure electrode drop and the voltage distribution across the electrode gap, nine 
pure tungsten wire probes (0. 76 mm diam), insulated electrically from the test section, 
are placed on both sides of the test section in the direction of the B field and in the plane 
of the current path. They are spaced equally in a 3. 8-centimeter electrode gap as 
shown in figure 6. 



Figure 6. - End view of test section showing electrodes and voltage probes. 
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Experimental Measurements 


: ' Arrm conditions. - The argon stagnation temperature was calculated calorimetric- 
ally. (Power into the gas is equal to the power into the arc minus the power into arc 
cooling water. ) The static temperature was determined by assuming an isentropic ex- 
pansion. The pressure was measured with a manometer. 

Sheath voltage. - Cathode sheath voltage measurements were made for voltages ap- 
plied across the electrodes and for voltages induced by flow through a magnetic field 
(i. e. , MHD generator). The 80-volt dc power supply used for the applied voltage is lim- 
ited to 50 amperes. The induced voltage is produced by a magnet with maximum mag- 
netic field strength of 0. 58 tesla. A typical set of voltage measurements is shown in fig- 
ure 7. The cathode sheath voltage drop was the same for the two voltage sources. 

Cathode spot temperature. - Electrode temperatures were measured by a two-color 
(wavelengths of 650 and 550 nm) pyrometer. The temperature is indicated on a meter 
calibrated in degrees with a linear scale. Also a millivolt output is used to record the 
data. 
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Figure S. - Test sections. 

The optical path for the pyrometer measurements included a quartz window and a 
mirror. Compared to a known temperature light source, this path introduced about a 
1 percent error in the temperature measurement. It was found that during a rim the 
quartz window did not discolor or lose its transparency. A cutaway sketch showing the 
pyrometer view and electrodes is shown in figure 8(a). 

Photographs. - Color movies (24 f/sec, f 22, 90 ND, Ektachrome 16 mm) were taken 
of the cathode spot. The current -voltage and electrode temperature data were recorded 
at the same time so that current, arc spot size and temperature could be correlated. A 
sketch showing the camera and electrodes is shown in figure 8(b). Enlargements of 
typical frames are shown in figure 9. 

Spectrographic study . - Spectrometer measurements were made for us by R. J. Sovie 
using a 1/2-meter Ebert monochrometer. A movable slit attached to a light-pipe scanned 
the gap between the electrodes. The lines examined were 434. 8 nanometers (ion line) and 
427. 2 nanometers (neutral line). The ion line appears at the cathode but not at the anode. 
This is shown in figure 10. The neutral line studied at the same time indicated a rela- 
tively uniform distribution in the electrode gap. 


Mirror 

Leads to 
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(a) Magnetic field applied, 20 amp. 



(b) No magnetic field, 40 amp. 

Figure 9. - Photographs of arc spots. Electroce size is 1.27x2.22 cm. Circles in 
center are support rod welds. 


18 




Upper electrode 


Lower electrode 


LZ 

(a) Upper electrode is cathode. 


Upper electrode 


Lower electrode 




Intensity ► 

(b) Lower electrode is cathode. 


Figure 10. - Argon ion line (434. 8 nm) intensity profile for discharge 
current of 30 amperes. Distance between electrodes, 1. 27 centi- 
meters. 


RESULTS AND DISCUSSION 
Comparison of Model and Experiment 

Presence of cathode spot and ionization region . - The sample photographs (movie 
frames) of cathode spots shown in figure 9 were for two currents. Note that the spot 
becomes brighter and larger at the higher current. (The bright circle in the center of 
the electrode is the support rod weld. ) 

From the ion line intensity profiles (shown in fig. 10), it is apparent that there is a 
region with large ion number density near the cathode, regardless of whether it is the 
upper or lower electrode. This region is identified with the ionization region of the model 
used in the analysis. 

Sheath voltage and c athode spot temperature . - Equation (32) provides a relation be- 
tween the spot temperature and the sheath voltage. Three parameters must be specified: 


TABLE I. - EXPERIMENTAL CONDITIONS IN TEST SECTION 


Run 

Pressure, 

p, 

N/cm^ 

Ion 

tempera- 

ture, 

T i- 

K 

Work 

function, 

eV 

Undisturbed 

cathode 

temperature, 

Tec, 

K 

Maximum 
ion cur- 
rent 
density, 
h, max’ 
A/cm 2 

1 

0.234 

990 

3.25 to 3.35 

2080 

499 

2 

.298 

540 

2.7 

1780 

858 

3 

.234 

1025 

3.2 to 3.6 

2080 

489 
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the argon pressure, argon temperature, and the work functions. These parameters are 
shown in table I for three experimental runs. The resulting theoretical values are com- 
pared with the experimental data in figure 11. The work function shown in table I was 
determined as the value which gave the best fit to the experimental data. Further dis- 
cussion on this point is found in the section entitled Work functio n. 

The curves seem to fit the data except when the sheath voltage is below 14 volts, 
which is near the excitation potential (11. 4 V). It is likely that the model describing the 



Figure 11. - Cathode sheath voltage as function of cathode spot temperature at 
two ion temperatures and three different values of work function at each ion 
temperature. 


processes in the ionization region is inaccurate for voltages near the excitation potential. 
Also, the spot radius becomes larger at high spot temperature and may not be small 
compared to the cathode dimensions. 

Current -vo ltage characteristics. - Using the parameters in table I, the current- 
voltage characteristics were calculated from equation (27), with 0 determined from 
equation (32). In figure 12 these characteristics are compared with the experimental 
data. A reasonable agreement is noted, giving added confidence to the values of work 
function selected from figure 11. 

Spot radius an d spot temperature. - The spot radius as a function of spot temperature 
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<b) Run 2. 


Figure 12. - Cathode sheath voltage as function of total current in 
argon for different values of cathode work function. 


calculated from equation (34) is compared in figure 13 with values measured from photo- 
graphs similar to those in figure 9. The radius is defined as 


R = 



(42) 


where the area is measured by using a planimeter. The agreement is again reasonable 
for the work functions determined from figure 11. Further evidence of agreement 
from the same experiment is the fact that the radius is quite sensitive to changes in spot 
temperature and work function. 

Work function . - The work function obtained from the best fit of the sheath voltage 
as a function of cathode spot temperature data (fig. 11) is plotted as a function of spot 
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Figure 13. - Cathode spot radius as function of spot temperatures at 
two different undisturbed cathode temperatures. 


temperature in figure 14. Also in figure 14 the work function is obtained from figure 13. 
The agreement between the two methods is reasonably good. The data used in the two 
methods was taken on different rims. 


From figure 14, the work function for the unactivated cathodes used in the experiment 

is 2.7. The Richardson equation constant in our calculations is set equal to 1. 2x10 
2 2 

amp/(m )(K ). However, other experimentors (e. g. , ref. 11) have found work functions 
for activated thoriated tungsten in vacuum to be 2. 6 eV and A . ^ equal to 3x10^ 
amp/(m 2 )(K 2 ) or 3. 2 eV using our assumed A rich of 1. 2xl0 6 amp/(m 2 )(K 2 ) at 1900 K. 

The difference of about 0. 5 eV in the work function mentioned above may be explained 
by the difference in cathode preparation (activation process) and cathode operating condi- 
tions. 

During the activation process (electrodes being heated to a high temperature in our 
case) oxides are reduced. The presence of oxygen on the cathode may explain our lower 
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Figure 14. - Work function determined from figures 11 and 13 as function of 
cathode spot temperature. 


work function once emission has been obtained. Many investigators have noted that both 
for tungsten cathodes in the presence of cesium and barium (ref. 10, pp. 56-58) and for 
copper electrodes in the presence of cesium (ref. 12) emission is improved (i.e. , the 
work function is lowered) in the presence of oxygen. Since "the mechanism of the emis- 
sion of the tungsten -thorium cathode is exactly the same as that of tungsten-cesium and 
the tungsten -barium cathodes. " (ref. 11, p. 55), it seems likely that oxygen has also 
lowered the work function of the unactivated thoriated tungsten electrodes tested herein, 
as well as permitting higher operating temperatures. It can also be seen from figure 14 
that at 2300 K the work function of our unactivated cathodes has increased to 3. 3 eV, in- 
dicating that the active material is being removed and that the work function is approach- 
ing that of pure tungsten, 4. 5 eV. 

Further evidence that the presence of oxygen lowered the work function of the cath- 
odes was gathered by carbonizing (as described in ref. 13) some cathodes and studying 
them. The work function at 2150 K increased from 3. 1 for the uncarbonized cathode to 
3. 5 eV for the carbonized cathode. Also, no spots were observed. The work function for 
lower cathode temperatures (i. e. , without a spot) was probably lowered by the activation 
process. However, at higher temperatures (i. e. , with a spot when using unactivated 
cathodes) the presence of oxygen for the unactivated cathode then provided the lower work 
function described previously. This explanation is consistent with the observations of 
Hermann and Wagner who state "that the work function of atomic film cathodes such as 
W-Ba, W-Ca, etc. can be both decreased and increased by the absorption of oxygen, and 
that the position of oxygen relative to the other absorbed material is decisive. " (ref. 11, 
p. 58). 
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Comparison with Previous Models 


The present model incorporates the properties of the cathode and the plasma in such 
a way that the cathode spot parameters can be deduced with no emperical constants re- 
quired. Therefore, the results from this model can be compared directly only with those 
of the Bade and Yos model. For the experimental conditions considered herein, the Bade 
and Yos model gives sheath voltages about 4 volts lower, current density about twice as 
large, a radius about 2/3 as large, heat transfer to the cathode about 30 percent less, 
and a spot temperature about 100° higher than the present model. The present model is 
an improvement of the Bade and Yos model without an undue increase in complexity since 
it has been shown that the present model agrees fairly well with the experiment. 


Implication for Generator Design 

It is desirable to maintain a low heat transfer to the cathode of an MHD generator. 
This is an energy loss and contributes to reducing the efficiency. In order to maintain 
a low heat transfer to the cathode the undisturbed cathode temperature should be kept 
high. These high temperatures also reduce the sheath voltage drop at the cathode, which 
is important in providing generators with high isentropic efficiency. 

The heat transfer also decreased with increasing pressure. However, the spot tem- 
perature became high and spot size became small. This may lead to a reduced cathode 
life because of erosion and depletion of thoria. 


CONCLUDING REMARKS 

The simple model used in this study agrees quite well with some experimental data 
for electrode phenomena. The results are qualitatively the same as that of Bade and Yos, 
but agreement with our experiment is improved. A major conclusion is that higher tem- 
perature electrodes will decrease the sheath voltage and, hence, the heat transfer losses 
to the cathode. This conclusion is important in determining the performance of MHD 
generators. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 23, 1969, 

129-02-08-05-22. 
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APPENDIX - SYMBOLS 


■^rich 

Richardson’s constant, 

e 

emissivity 

1.2X10 6 A/(m 2 )(K 2 ) 

9 

nondimensional temperature 

e 

electrical charge, 1. 6x10” a C 

X 

thickness of ionization region 

I 

current 

o 

Stefan-Boltzmann constant, 

i 

current density 


5. 75X10 -12 W/(cm 2 )(K 4 ) 

K 

thermal conductivity 

$ 

parameter defined in eq. (25) 

k 

Boltzman constant, 
1. 38X10" 23 J/K 

<P 

potential (ionization, excitation) 
or work function 

m 

particle mass 

* 

parameter defined in eq. (29) 

N 

number density 

Subscripts: 

P 

pressure 

e 

electron 

P 

parameter defined in eq. (15) 

i 

ion 

Q 

heat transfer 

o 

neutral 

R 

spot radius 

ref 

reference 

r 

ratio of electron to ion temper- 

s 

spot 


ature 

sh 

sheath 

T 

temperature 

w 

work function 

V 

voltage 

1 

excited state 

v i 

velocity 

oo 

at a large distance from spot 
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